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ABSTRACT. The Cauchy problem for the Korteweg de Vries (KdV) equation with small 
dispersion of order e 2 , e <C 1, is characterized by the appearance of a zone of rapid modu- 
lated oscillations. These oscillations are approximately described by the elliptic solution of 
KdV where the amplitude, wave-number and frequency are not constant but evolve accord- 
ing to the Whitham equations. Whereas the difference between the KdV and the asymptotic 
solution decreases as e in the interior of the Whitham oscillatory zone, it is known to be 
only of order e 1 ' 3 near the leading edge of this zone. To obtain a more accurate descrip- 
tion near the leading edge of the oscillatory zone we present a multiscale expansion of the 
solution of KdV in terms of the Hastings-McLeod solution of the Painleve-II equation. We 
show numerically that the resulting multiscale solution approximates the KdV solution, in 
the small dispersion limit, to the order e 2 / 3 . 



1. Introduction 

The mathematically rigorous study of the small dispersion limit of the Korteweg de 
Vries (KdV) equation 

(1.1) u t + 6uu x + e 2 u xxx = 0, £<Cl, 

with smooth initial data uq(x) was initiated in the works of Lax-Levermore [28 1, which 
stimulated intense research both numerically and analytically on the problem. The solution 




Figure 1 . Numerical solution of the KdV equation for the initial data 
uq(x) = — sech 2 x and e = 0.1. 
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of the Cauchy problem of the KdV equation in the small dispersion limit is characterized 
by the appearance of a zone of rapid oscillations of frequency of order 1/e, see for instance 
Fig.Q] 

These oscillations are formed in the strong nonlinear regime and they have been ana- 
lytically described in terms of elliptic functions, and in the general case in terms of theta 
functions in [34], [8 |; the evolution in time of the oscillations was studied in [ 30 1 . These 
results give a good asymptotic description of the oscillations only near the center of the 
oscillatory zone (see Fig. 13. In lfl9l we have studied numerically the small dispersion 
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Figure 2. We plot for u = -sech 2 a;, t = 0.4 and e = 10~ 2 from 
top to bottom: 1) the numerical solution of KdV; 2) the asymptotic for- 
mula ( 12. 3t in terms of elliptic functions and the Hopf solution; 3) the 
multiscale solution where the envelope of the oscillations is given by a 
solution to the Painleve-II equation. 

limit of the KdV equation for the concrete example of initial data of the form 
(1.2) uq(x) = — sech 2 ir. 

We have compared the asymptotic description given in the works l28l [34l [8l with the 
numerical KdV solution. In lfl9l we have shown numerically that the difference between 
the KdV solution and the elliptic asymptotic solution at the center of the oscillatory zone 
scales like e while this fails to be true at the boundary of the oscillatory zone. This fact 
was also observed for the Benjamin-Ono in |24l . In particular at the left boundary, where 
the amplitude of the oscillations tend to zero, the difference between the KdV solution 
and the elliptic asymptotic solution scales like ea. In this manuscript we show that the 
Painleve-II equation describes the envelope of the oscillations at the leading edge where 
the oscillations tend to zero. More precisely 

u(x, t, e) = u(t) + e^a cos f t\ + 0(e 2 / 3 ) 
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Figure 3 . In (a) the difference between the upper two plots of Fig. |2]is 
shown. The Whitham zone is shown in blue. In (b) one can see the same 
situation as in (a) except for the region close to the leading edge of the 
Whitham zone where the difference between the KdV and the multiscale 
solution is shown in red. 



where u(t) — u(x~ (t), t) solve the Hopf equation at the boundary x~ (t) of the oscillatory 
zone, ip is a phase determined in ( I3.291 l and a is, up to shift and rescalings, the Hastings- 
McLeod solution of Painleve-II equation 

A zz = zA + 2 A 3 

with boundary conditions 

lim A{z) = yJ-z/2, lim A(z) = Ai(z) 

z — > — oo z — >+oo 

where Ai(z) is the Airy function. 

Painleve equations appear in many branches of mathematics (for a review see [6|). For 
example in the study of self-similar solutions of integrable equations, in the study of the 
Hele-Shaw flow near singularities 04), or in double-scaling limits in random matrix mod- 
els (see e.g. (T), Ifl5l .ll2l.ll5l). In this work, following [27|, we perform a double-scaling 
limit of the KdV equation to derive the asymptotic description of the leading edge oscil- 
lations which are formed in the KdV small dispersion limit. We show that the envelope 
of the oscillations is determined by the Hastings-McLeod [22] solution of the Painleve- 
II equation. Then we compare numerically for the initial data uq(x) = — sech 2 x at the 
leading edge of the oscillatory front, the KdV solution with the derived multiscale solu- 
tion and show that the difference between the two solutions scales like ea . We identify a 
neighborhood of the leading edge of the Whitham zone where the multiscale solution gives 
a better asymptotic description than the asymptotic solution based on the elliptic and the 
Hopf solution. This allows to patch different asymptotic descriptions to provide a more 
satisfactory treatment of the small dispersion limit of KdV as shown in Fig. [3] 
Our analytical investigation of the multiscale expansion of the KdV solution requires the 
following assumptions on the initial data: 

• uq (x) is negative with a single minimum; we chose the minimum to be at x = 
and normalized to -1, namely uo(0) = —1; 

• the function /_ (u) which is the inverse of the monotone decreasing part of the 
initial data uo(x) is such that f'l'(u) < 0; 

• uq(x)(1 + x 2 )dx < oo. 
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The last condition guarantees that the solution of the Cauchy problem for KdV exists for 
all times t > 0. 

This manuscript is organized as follows. In section 2 we review the small dispersion 
limit of the KdV equation and study the small amplitude limit of the oscillations. In section 
3 we perform a multiscale expansion of the KdV solution when the oscillations tend to zero, 
and we show that the envelope of the oscillations is given by a solution of the Painleve- 
II equation. In section 4 we numerically compare the KdV solution with the multiscale 
solution obtained in section 2. We show that the difference between the two solutions 

2 

scales as which is in accordance with our analytical result. We identify a zone near the 
leading edge where the multiscale solution provides a better description than the elliptic 
and the Hopf solution and patch the solutions. In section 5 we summarize the results and 
add some concluding remarks on future directions of research. 

2. Asymptotic solution of KdV in the small dispersion limit 

The solution of the Cauchy problem for the KdV equation is characterized by the ap- 
pearance of a zone of fast oscillations of wave-length of order e, see e.g. Fig. Q] These 
oscillations were called by Gurevich and Pitaevski dispersive shock waves [21 1. Following 
the work of [28], ll34l and [8|, the description of the small dispersion limit of the KdV 
equation is the following: 

1) for < t < t c , where t c is a critical time, the solution u(x, t, e) of the KdV Cauchy 
problem is approximated, for small e, by u(x, t) which solves the Hopf equation 

(2.1) ut + 6uu x = 0. 

Here t c is the time when the first point of gradient catastrophe appears in the solution 

(2.2) «(x,t)=uo(0, x = Gtuo(0+t, 

of the Hopf equation. From the above, the time t c of gradient catastrophe can be evaluated 
from the relation 

^ _ 1 

max 5eR [-6-u^(C)]' 

2) After the time of gradient catastrophe, the solution of the KdV equation is characterized 
by the appearance of an interval of rapid modulated oscillations. According to the Lax- 
Levermore theory, the interval [x~ (t), x + (t)] of the oscillatory zone is independent of 
e. Here x~(t) and x + (t) are determined from the initial data and satisfy the condition 
x~{t c ) — x + (t c ) — x c where x c is the x-coordinate of the point of gradient catastrophe 
of the Hopf solution. Outside the interval [x~(t), x + (t)] the leading order asymptotics 
of u(x, t, e) as e — > is described by the solution of the Hopf equation (12.21 i. Inside the 
interval [x~(t), x + (t)] the solution u(x,t,e) is approximately described, for small e, by 
the eUiptic solution of KdV ED , J28) , E) , El , 

(2.3) u(x, t, e) ~ +fa +(3 2 + (33 + 2a + 2e 2 -^ log B(Sl(x, t);T) 
where 

(2.4) n = V 2eK{3) 3 [X ~ 2t{Pl + ^+^~^ 
and 

(2.5) a = -Pi + (Pi-P 3 )——, T = i——, s =- — 

K(s) K{s) Pi- Ps 

with K{s) and E(s) the complete elliptic integrals of the first and second kind, K'{s) — 
K{s/1 — s 2 ); 9 is the Jacobi elliptic theta function defined by the Fourier series 



(2.6) 6{z-T) = Y, 



^7rin T-\-27rinz 
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For constant values of the fa the formula d2.31 > is an exact solution of KdV well known in 
the theory of finite gap integration l23l . iflOl l. However in the description of the leading 
order asymptotics of u(x, t, e) as e — * 0, the quantities fa depend on x and t and evolve 
according to the Whitham equations 



(2.7) -f3 i + Vi —(3 i = ^ i= 1,2,3, 
where the speeds Vi are given by the formula 

(2.8) Vi = 4 n ^- /3fc) + 2(ft + fa + fa), 

with a as in (12.5l l. The formula for q in the phase 51 in (12.4l i that we are giving below was 
introduced in [ 19] and looks different but is equivalent to the one in j8) 

i r' r 1 , i-O^i^jh + MA) + 



as, ^.A.«- 5 5 K y. i y. i ** 

where f-(y) is the inverse function of the decreasing part of the initial data no. The above 
formula for q(fa, fa, fa) is valid as long as fa > fa > fa > — 1- When pa reaches the 
minimum value — 1 and passes over the negative hump, it is necessary to take into account 
also the increasing part of the initial data f+(u) in formula ( 12.91 ). We denote by T this 
time. For t > T > t c we introduce the variable X$ defined by ^(X^) = fa which is still 
monotonous. For values of X$ beyond the hump, namely X% > 0, we have to substitute 
i by the formula 




(2.10) q(fa,fa,fa)-- 

The function q — q{fa,fa, fa) is symmetric with respect to fa, fa and fa, and satisfies a 
linear over-determined system of Euler-Poisson-Darboux type. It has been introduced in 
the work of Fei-Ran Tian 1 30 1 . The Whitham equations d2.7l) can be integrated through the 
so called hodograph transform, which generalizes the method of characteristics, and which 
gives the solution in the implicit form ll33l 

(2.11) x — v{t + u>i, i — 1,2,3, 

where the Vi are defined in (12.8b and the Wi = Wi (p\ , fa , fa ) are obtained from an algebro- 
geometric procedure [26 1 by the formula | 

(2. 12) Wi = l(vi-2Y0k I - - 'I- ' = 1-2. 3. 




with q defined in i2.9\ or (12.10b . Formula ( 12.1 U solves the initial value problem for the 
Whitham equations ( 12.7b with boundary conditions: 

a) leading edge: 

fa = the Hopf solution ( 12. lb 

(2.13) 

P2 — P3, 

b) trailing edge: 

fa=fa 

(2 14) 

fa = the Hopf solution dm. 

In lfl9l we have solved numerically the initial value problem for the Whitham equations. 
In this way we could perform a numerical comparison between the KdV small dispersion 
solution and the asymptotic formula d2.3b (see Fig. [3). While in the interior of the oscil- 
latory zone the error scales numerically like e, at the left boundary of the oscillatory zone 
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the error scales numerically like es. To derive a more satisfactory asymptotic approxi- 
mation of the KdV small dispersion limit in the vicinity of this point, we perform in the 
next section a double scaling expansion of the KdV equation, following the double scaling 
limits appearing in random matrix theory. Before doing this analysis, we study the elliptic 
solution (12.3b in the limit when the oscillations go to zero. 

2.1. Small amplitude limit of the elliptic solution. We study the elliptic solution (12.31 i 
near the leading edge, namely when oscillations go to zero. To avoid degeneracies, we 
rewrite the system ( 12.1 lb in the equivalent form 



(vit + wi -x)(a + (3 1 ) = 

1>2< + W2 — X = 
1 



\(v 2 - v 3 )t + w 2 - w 3 ] = 0. 



(2.15) 

(ft - ft 

and perform the limit 5^0 where 

{h = v + 5, /3 3 = v — 5, ft = u. 

To simplify our calculation we restrict ourselves to the case t c < t < T. The following 
limit holds: 

K(s) v — pi i(v — pi) z 

such that 
(2.17) 



A(u-v)' 

Furthermore the following identities hold 

(2.18) f-{u) =[2(u - v)d u q(u,v,v) + q(u,v,v)} 

(2.19) u) =d v q(u, v, v) + d u q(u, v, v) 
where 

(2.20) ^(v,u) = — ] = , =* — - = —, / dfi-~ 



2y/2J-l V 1 - I 1 %VV - U J u yjv - fl' 

Substituting (12. 17b (12. 1 8b and ( 12.191 ) into ( 12.151 ) we arrive at the system 
(2.21) 

2 (x-6tu- f-(u)-2(u-v)(6t + $(v;u)) 4 
x = btu + j-{u) — — ^5 h O(0 ), 

H(v — uy 

x = 6tu + /_ (u) + 2(v - u)[6t + $(t>, u)] + 5[6t + $(t>, u) + (v - u)d v $(v, u)] 



5 2 3 

[6t - 2(u - v) 2 d vvv q(u, v, v) + 4(u - v)d vv q(u, v, v) + -d v q(u, v, v, )] + 0(6 3 ) 



= 6< + u) + (v - u)d v $(v, u) + 0(5). 

From the above we deduce that, in the limit 5 — ► 0, the hodograph transform ( 12.15b reduces 
to the form (see ll30llfT8l ) 

' 6ut + f-(u)-x = 
(2.22) I $(v,u) + 6t = Q 

^ d v $(v,u) = 0. 

The above system enables one to determine x, u and v as a function of time. This time 
dependence will be denoted x = x~(t), u = u(t) and v — v(t). We are interested in 
studying the behavior of the elliptic solution (12.3b near the leading edge, namely when 



MULTISCALE EXPANSION OF KDV 7 

x — x~ (t) is small and x > x~ (t). For this purpose we introduce two unknown functions 
of x and t, 

8 = 8(x-x-(t)), A = A(x - x~(t)) 

which tend to zero as x — ► x~(t). We are going to derive the dependence of A as a 
function of x — x~ (t). Let us fix 

(2.23) ft = v + S, p s = v-8, 5^0 ft = u + A, A -> 0. 

Using the first equation of (12. 2U we obtain 

= x - 6tu - /_(ft) + ^(-^-/-^)- 2 ^-^ + ^^)) + 0(O . 

8(ft - u) 2 

Expanding the above expression near ft (x, i) = u(t) + A(x, t), using the identity 

£«^>- '^--;r' 

and ( 12.221 1 we arrive at the expression 

(2.25) 0~x-x-(0-(6t + /iM)A+ * (x-x-(t)) ; 

so that 

x — x~ (t) 



(2.26) A_ , 

U + f_{u) 

Using the second equation in ( 12.21b we arrive at 

(2.27) x-x~(t)~8 2 c 

where 
(2.28) 

3 

c = [6i - 2(« - v) d vvv q(u, v, v) + 4(u - v)d vv q(u, v. v) + -d v q(u 7 v, v, )]/4(u - u) 

= — 

Therefore 

Theorem 2.1. The elliptic solution ( 12.31 ) in the limit (12.23\) takes the form 



, _x cos I 4tt- 



(2.29) M (a;,t,e)~ M (t)+ ^ + 25 cos 2tt— + J'- , 

where the phase Q~ takes the form 

(2.30) 2-kQT =0 o + 02 
with 



(2.31) o (t) =-16 / (u(t) -«(r))»dr, ^a(x,i) = 2V«(*) - «(i)(a? -a;-(t)), 



a«t/ u(t), ant/ x (t) soZve the system K2.22\ . 

Proof. We first prove the relation (12.30b . Using the expansion 
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and ( 12.18b we obtain the following limit for the phase £1 in ( 12.41 ) 

27^,3=^ = 2y/p 1 -v{l - ^——^ )[x - Gtpx - /_(&) 

+ 2(A - v)(2t + d 01 q(p u v,v) - § -dU(Pi,v, v)} + 0(5 4 ). 

Using the identity 

(I2.22l i and (12.19b we can rewrite the above in the form 

a2 (a x 3(2t + a^ 1 g(/3 1 ,t;,i;)) 3 6£ + /:( U ) A ^, A ^ 
»„?(ft,«,«) = ^T^) + 4 {v-uY ] 

so that the phase il takes the form 

3 5 2 x - x~(t) 

(2.32) 2^\ P2 3=v±s ~ 4(A - w)*(2* + d^q^v, «)) - V . 

° (« — U) 2 

We define 

J7o(A,«) :=VA - «[2(/3i - u)t + (A - v)d 01 q((3x, v, v)] 

(2.33) =2/ ^/ft -A[$(A,#t) + 6i]dA, 

so that, expanding near j3\ = u + A, by ( 12.22b 

(2.34) Vo(lh,v) = m (u, v) + 2 (ar, t) + 0(A 2 ), 

where <f>2(x, t) is defined in ( 12.311 ). To show that rjo(u, v) defined in (12.33b coincides with 
4>o defined in ( 12.31b . we differentiate ( 12.33b with respect to time, 

d 



2— / Vu- A[$(A, u) + 6t]d\ = -16(u - v)* 

J v 

(u — n\ 

where we have used the identity ( 12.24b and dtuit) = 12 . Integrating the r.h.s 

6t + fL(u) 

of the above expression with respect to t from t c to t we obtain the formula ( 12.31b . 
Using ( 12.34b and ( 12.261 ) we rewrite the phase ( 12.321 ) in the form 

(2.35) 2^| fJl=u+A - 0o + 02 - * 2( * " X }P ] + 0(A 2 ) 

P 2 ,3=V±8 8{U — V)2 

where 0o and 02 are as in ( 12.311 ). Neglecting the higher order terms in 5 and A of the 
above expansion one obtains (12.30b . 

Now we are ready to expand the theta-function expression in the limit of small ampli- 
tudes. Using ( 12.171 ) and 

e" T - sr^Cl - — ) + 4t-Aj + 0(6 5 ,5 S A), 
8{u — v) u — v 128 [u — w) J 

one derives the small amplitude limit of the Jacobi 6>-function ( 12.6b 

00; r) = 1 + - - - (1 - -^—) cos(27rz) + 0(5 4 ). 
4(m — v) u — v 

Substituting the above expansion in (12.3b one obtains 



x — x~(t) ( Q 
u(x, t, e) ~ u(t) + — -jA- + 25 (x, t) cos 2?r I + 5(x, t) 



cos 47r — 1 



6t + f'_(u) v ' ' V e / ' 2(u(t) -«(*)) ' 
which coincides with ( 12.29b . □ 
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3. Painleve equations at the leading edge 

In this section we present a multiscale description of the oscillatory behavior of a solu- 
tion to the KdV equation in the small dispersion limit close to the leading edge x~ (t) where 
/?2 = 03 = v and (3i = u. We are interested in the double scaling limit to the solution of 
the KdV equation (11. Il l as e — > and x — ► x~ (t) in such a way that x — x~ (t) oc e 2 / 3 . We 
introduce the rescaled coordinate y near the leading edge, 

(3.1) y = e- 2 / 3 (x-x-(t)), 
which transforms the KdV equation (11. j} , to the form 

(3.2) e3U t + e^Uyyy + (6u — a; t )% = 0, 

where = ~^ x ^ (*)■ The substitution ( 13. U has the effect that the linear term of (13.2l i is 
just the Airy equation u t + u yyy = which has oscillatory solutions. 

It is known [9|,|27| that the corrections to the Hopf solution near the leading edge are 
of the order e 1 / 3 . We thus make the ansatz 

(3.3) u(y, t, e) = U Q + e 1/3 C/i + e 2/3 C/ 2 +eU 3 + ..., 

where Uo = u(t) is the solution at the leading edge. We assume that Uk>i contains 
oscillatory terms with oscillations of the order 1/e. In particular 



(3.4) U 1 = a(y,t) cos 



where 

(3.5) i>( y , t) = My, t) + e 1/3 Mv, t) + e 2/ V 2 (y, t) + e^ 3 (y, t) 

Similarly we put 

(3.6) U 2 = b 1 (y,t) + b 2 (y,t) cos' 2r[ " J) 
and 

(3.7) U 3 = co(y, t) + c 2 (y,t) sin I 1 + c 3 (y,t) cos I 

Terms proportional to sin(ip/e) can be absorbed by a redefinition of tp. Since we impose 
no further restrictions on ip here, such terms are therefore omitted in all orders. We only 
consider terms proportional to cos(ip/e) in order e 1 / 3 and the necessary terms in higher 
order to compensate the terms due to the nonlinearities in (13.21 i. 

If we enter equation ( 11. U with this ansatz, we immediately obtain from the term of 
order e° that ipo !V — fa y — 0. From the term of order we get 

(3.8) i>l y -(6U -xi)fa y -faj =0. 
In order e 2 / 3 we obtain the following equations 

<3 - 9) 

(3.10) fayiWty - QU + x^)- fa it - 

(3.11) ^-[a 2 (3^l y -6U Q + x;)} = 0. 
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In order e we get 

(3.12) U ,t + (6U Q - Xt)b hy + 3aa y = 

(3.13) i> 2 , y ^ {a 2 ip 3 , y ) = 

(3.14) 2a 2 ih lV = 

(3-15) 1$, v OvV + a ( 2b l^ly + l^2,yMt) + \a 3 = 



•-> aa v 
(3.16) c 2 = — -f- 

V2,y 



3a 



(3.17) c, , 

16^2,2/ 

A solution to dllOl . dTTTT i and (13141 is 

(3.18) 3^,=6[/ -x- + ^, Vi,t = 0, 

where C(t) is a free function of i. Note that 

(3.19) U (t) = u(t), Xt = I2v{t) - 6u(i), 

where u(t) and v(t) are defined in ( 12.22b . The above equations imply that 

Mv,t) = 2yV(u-v) + J ^Ld y + h(t), 

where h(t) is free function of t. Comparing the above relation with the formula of 4>2 in 
(12. 3U of the phase in the small amplitude expansion we can conclude that C(t) = and 
h(t) = so that 



(3.20) MV, t) = 2yy/u(t)-v{t). 
From d3~8b . ( f3TT9b and ( f3~20l ) we derive that 

V>cu = — 16(u — v)*, 

namely 

(3.21) Mt) = -16 / (u(r) - u(r))ldr. 



From (I3TT21 we find 

a?_ yUo_ 



&1 = - 0T2 ^7l2 ^ fc W 



(3.22) 

a y , . 

= ~8[u(i) + 6t+7T(u) + W ' 

where fc(t) is a free function of t. It will be fixed by matching it with the elliptic solution 
in the Whitham zone. Substituting ( I3.22l i and ( 13.20b into (13.15b we arrive at the equation 

,s 2 . . / 12k(t)(u(t) - v(t))\ a 3 

(3.23) 4(«(t) - «(t))a w - -«t(t)a U = y. 

Making the substitution 

(3.24) A = 6*a/(4 J/» (u _ t , ) x/- )j z= _ ) ^ (y _ , Q ) 

with ?/o = 12fc(t)(u(t) — v(t))/vt, we arrive at the equation 

(3.25) ^ = + 

which is a special case of the Painleve-II equation A zz = zA + 2A 3 — 7, with 7 = 0. 
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Since we are only interested in terms up to order e 1 ' 3 in u(x, t, e), the terms b\, 62, cq, 
c-2 and C3 are not important for us. However, we had to go to order e to determine ip^ which 
will contribute to the e 1 / 3 terms in u. 

To sum up we get for u(x, t, e) 
(3.26) 

"a 2 (cos(2V>/e) - 1) 



u(x, t, e) = u(t)+e 1/3 a cos f ^ J +e 2/3 
where 



8(«(i) -«(*)) 



•*(*) 



2/ 

6t + /i(«) 



0(e) 



V»(y,i) = -16 / (u(r) - v(t))Ut + e*h + e%[2yy/u(t) - v(t)} + etp 3 {t) + 00), 



and h is an integration constant. There are free functions of t in the integration of the multi- 
scale equations, namely the functions k(t) and ipz(t) and the constant h. Moreover, the 
solution of the Painleve II equations needs to be fixed. We fix the constants by comparing 
the multiscale solution ( 13.281 ) with the elliptic solution ( 12.29b at the border of the Whitham 
zone. Indeed comparing ( 13.26b . ( 12.29b and ( I2.351 l we obtain 



(3.27) 



1 ! 2 

S = 2 e3 °' A = 0(e=>' 



and k(t) — 0, ip3(t) = and h = 0. Therefore the multiscale solution takes the form 
(3.28) 



u(x, t, e) = u(t)+e 1 ^ 3 a(y, t) cos 



a 2 (cos(2V>/e) - 1) 



y 



8(u{t)-v{t)) 6t + fL(u) 
where y = e~ 2 / 3 (x — x~{t)), a(y,t) satisfies (13.231 1 and 



ip(y,t) = —16 / (u(t) — v 



(3.29) 



For the numerical comparison in the following section, we consider terms up to order e 1 / 3 
in u in ( 13.28b 



(3.30) 



u(x, t, e) = u(t) + e^aiy, t) cos ft^Tj + 0(e 2 / 3 ), 



where ip(y,t) is as given above. For fixing the particular solution of the Painleve-II 
equation (13.251 l the following considerations are needed. For large x < x~(t), the so- 
lution of KdV is essentially approximated by the solution of the Hopf equation, and the 
term of order e 1 / 3 has to be negligible in (13.30b . namely a(y) ~ for large negative 

y = (x — x~(i))e~ 2 / 3 . For 2 < x~ (t) 

1/3 



( - , , Vt , v \ ' (x-x-(t))e- 2 ^>0 



(> 



because v t = - — 

(u - v)d vv <s>(v;u) 

x <C x~(t) we conclude that 
(3.31) 



< since d vv Q(v; u) < and u > v. Therefore for 



lim A(z) = 0. 

z — >+oo 



For x > x (t), from the small amplitude limit of the elliptic solution of the KdV equation 
we obtain combining ( 12.27b . ( 12.28b and ( 13. 27b 



2{x - x-(t)) 
(u - v)d vv $(v;u) 
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which, in the limit e — * or y — > +00 gives 



lim a(y) 



I!/ 



y^+oc y (u — v)d vv Q(v; it) 

or equivalently, by (13.241 1. 



(3.32) lim A(z) = yJ-z/2. 



z — > — oc 



The existence and uniqueness of the solution of (13.25b satisfying (13.3 It and ( 13.32b was first 
established by Hastings and Mcleod [22| (see also later works of [25 1 and |6|). It is also 
worth noticing that the asymptotics of A(z) at z — > +00 can be specified as 

(3.33) lim A(z) = M(z) 

z — >+oo 

where Ai(z) is the Airy function. Moreover the asymptotic condition ( 13.331 ) characterizes 
the solution A(z) uniquely, so that ( 13.32b and ( 13.331 ) constitutes an example of the so called 
connection formula for the Painleve equations (see e.g. [ 13|, |7|). 



4. Comparison of the multiscale expansion and the asymptotic solution 

TO THE SMALL DISPERSION KdV 

The numerical evaluation of the asymptotic solution based on the Hopf and the elliptic 
solution is described in [19]. To evaluate the multiscale solution ( 13.281 ), one needs in 
addition to the quantities computed there the Hastings-McLeod solution to the Painleve- 
II equation. This solution was calculated numerically by Tracy and Widom ll32ll with 
standard solvers for ordinary differential equation and by Prahofer and Spohn ll29l [36l 
with in principle arbitrary precision with a Taylor series approach. The general family of 
solutions of Painleve-II such that lim z ^ +QO A(z) = aAi(z) with a positive constant was 
studied numerically in OTI and analytically in [7[. Solutions to Painleve-II in the complex 
plane were studied analytically and numerically by Fokas and Tanveer in ipHl . We use 
here an approach based on spectral methods which is described briefly in the appendix. 
This approach is both efficient and of high precision and can directly be combined with the 
numerics of [ 19]. 

Times t ^> t c . Close to breakup the multiscale expansion is expected to be inefficient since 
it is best near the leading edge, and since at breakup both the leading and the trailing edge 
coincide. We will discuss this solution close to breakup below, but first we will study it 
for time t = 0.4 ^> t c — 0.216 .... In Fig. H]one can see that the multiscale solution 
gives an excellent approximation of the KdV solution for x < x~ (0.4) = —3.2297 and in 
the Whitham zone close to x~ . For larger values of x, the solutions are out of phase and 
the values of the multiscale solution are shifted towards positive values. The difference 
of the two solutions is shown in Fig. [5] From this figure it is even more obvious that the 
multiscale solution is a valid approximation in the Whitham zone near the leading edge, 
but the difference increases rapidly for |x| 3> x~ . 

e dependence. In [|T9| it was shown that the asymptotic description becomes more accurate 
with decreasing e. The same is true for the multiscale solution as can be seen in Fig. [6] 
The zone, where the multiscale solution gives a better approximation than the asymptotic 
elliptic solution, shrinks with e. For x ^> x~(t), the multiscale solution is always only 
a poor approximation to the KdV solution. The maximal difference A max of the KdV 
solution and the multiscale solution near this edge decreases roughly as e 2 / 3 . More pre- 
cisely the error can be fitted with a straight line by a standard linear regression analysis, 
— log 10 A max — — alog 10 e + b with a = 0.63, b = 0.41. The correlation coefficient is 
r = 0.999, the standard error is a a = 0.02. 
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Figure 4. The blue line is the solution of the KdV equation for the 
initial data uq(x) — — sech 2 x and e = 10~ 2 for t = 0.4, and the green 
line is the corresponding multiscale solution given by formula ( 13.28b . 



0.1 p 
0.08 - 




Figure 5. The difference of the KdV and the multiscale solution for 
the initial data uq(x) = sech 2 x and e = 1CP 2 for t — 0.4. The curve is 
plotted in green in the Whitham zone. 

Comparison and matching with the asymptotic solution. The aim of this paper is to im- 
prove the asymptotic description of the small dispersion limit of KdV near the leading 
edge. In Fig. [7] it can be seen that the multiscale solution will indeed be a much better 
approximation near this edge. Near the leading edge, the multiscale solution provides a 
superior description of the KdV solution, whereas the elliptic asymptotic solution is much 
better for x ^> x~(t) in the Whitham zone. In fact it is possible to identify a zone where 
the multiscale solution is more satisfactory than the asymptotic solution. Due to the strong 
oscillations of the solutions, there is a certain ambiguity in the definition of this zone. We 
define the limits of the zone as the last intersection (or where the solutions come closest) 
on which the other solution has an error with larger oscillations. In this zone it is possi- 
ble to replace the asymptotic solution by the multiscale solution. The result of this patch 
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Figure 6. Difference of the KdV and the multiscale solution in order 
e 1 / 3 for the initial data uq(x) — — sech 2 x and several values of e for 
t = 0.4. The curves are plotted in green in the Whitham zone. 
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Figure 7 . The difference of the KdV and the multiscale solution (blue) 
and the difference of the KdV and the asymptotic solution (green) for the 
initial data uo(x) = — sech 2 a; and e = 10 -2 for t — 0.4. 

work approach is shown in Fig. [8] It can be seen that the resulting amended asymptotic 
description has an accuracy near the leading edge of the same order as in the interior of the 
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Figure 8 . Difference of the KdV and the multiscale solution (blue) and 
the KdV and the asymptotic solution (green) for the initial data uq (x) = 
— sech 2 a: at t = 0.4 for two values of e. 

Whitham zone. The maximal difference between the KdV and the asymptotic solution still 
occurs near the leading edge. 

As already mentioned, the zone where the multiscale solution provides a better approx- 
imation to the KdV solution than the asymptotic solution, shrinks with e as can be inferred 
from Fig. [9] The width of this zone decreases roughly as e 2 / 3 which shows the self consis- 
tency of the used rescaling of the spatial coordinate near the leading edge. More precisely, 
we find a scaling e a with a = 0.66, correlation coefficient r = 0.9996 and standard error 
a a = 0.015. It can be seen that the zone is not symmetric around the leading edge, it 
extends much further into the Hopf region than in the Whitham zone. This is due to the 
fact that the multiscale solution is quickly out of phase with the rapid oscillations in the 
Whitham zone, and that the Hopf solution does not have oscillations. 
Breakup time. In |fl9l it was shown that the elliptic asymptotic solution is worst near the 
breakup of the Hopf solution. The multiscale expansion obtained in the previous section 
is not defined for times before t c , and it will be worst there, since it can be understood 
as an expansion around the leading edge of the Whitham zone. At breakup, however, 
leading and trailing edge coincide. Thus the approximation is rather crude there, but it 
increases in quality with time as can be seen in Fig.[l0] It is, however, interesting to study 
at which times the multiscale solution starts to give a better asymptotic description than 
other approaches. Dubrovin conjectured ifTTI that the asymptotic behavior of the KdV 
solution close to the breakup of the corresponding Hopf solution is given by a particular 
solution to the second equation in the Painleve-I hierarchy. In [20| we provided strong 
numerical evidence for the validity of this conjecture. The natural question is whether 
Painleve-I2 description near the critical point provides a satisfactory asymptotic solution 
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Figure 9. Boundary values of the zone where the multiscale solution 
provides a better approximation to the KdV solution than the asymptotic 
solution. The x-values of the boundary of this zone (normalized by x~) 
for t = 0.4 are shown for several values of e. 

for KdV till times where the the multiscale solution studied in the present paper provides 
a valid description near the leading edge. A comparison of Fig. [10] with a similar figure in 
ll20l shows that this is qualitatively the case. In Fig.[l0]the multiscale solution is shown for 
x < x + (t). Near breakup the approximation is only acceptable close to the breakup point. 
For larger times, more and more oscillations are satisfactorily reproduced by the multiscale 
solution. But the multiscale solution will only be a better approximation of the oscillations 
than the Hopf solution for times t 3> t c . 

5. Outlook 

In the present work we have considered a multiscale solution to the KdV equation in 
the small dispersion limit close to the leading edge of the oscillatory zone. We studied the 
solution up to order e 1 ' 3 and show that the envelope of the oscillations is described by the 
Hasting-McLeod solution of the Painleve-II. The validity of the approach in the considered 
limit was shown numerically. The double scaling expansion of the KdV solution in the 
small dispersion limit will be investigated with the Riemann-Hilbert approach and steepest 
descent method for oscillatory Riemann-Hilbert problem as done in [8]. The Riemann- 
Hilbert approach seems so far the only analytical tool to study the double-scaling expansion 
to the Cauchy problem of the KdV equation. This project will be the subject of our future 
research. 

As can be seen from Figure [5J the asymptotic solution of the KdV equation does not 
give a satisfactory description of the KdV small dispersion limit also at the trailing edge of 
the oscillatory zone. This problem will be investigated in a subsequent publication, too. A 
similar problem was tackled in the context of matrix-models in PI . 

Appendix A. Numerical solution of the Painleve-II equation 

We are interested in the numerical computation of the Hastings-McLeod solution to the 
Painleve-II equation 

(1.1) P H A := A zz -zA- 2A 3 = 
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t=0.218 



t=0.220 



t=0.222 



=> -0.5 




= -0.5 




-1.7 -1.65 -1.6 -1.55 -1.5 -1.7 -1.65 -1.6 -1.55 -1.5 -1.7 -1.65 -1.6 -1.55 -1.5 

X = m2A t=6 < 226 t=tf.228 

0i 1 1 



a -0.5 




=> -0.5 




=> -0.5 




-1.7 -1.65 -1.6 -1.55 -1.5 
1=0*230 



=> -0.5 




-1.7 -1 .65 -1 .6 -1 .55 -1 .5 -1 .7 -1 .65 -1.6 -1 .55 -1 .5 

x x 

t=0.232 t=0.234 

0i 1 Or 



=> -0.5 





-1.7 -1.65 -1.6 -1.55 -1.5 



-1.7 -1.65 -1.6 -1.55 -1.5 



-1.7 -1.65 -1.6 -1.55 -1.5 



Figure 10. The blue line is the solution of the KdV equation for the 
initial data uq(x) = — sech 2 x and e = 10~ 2 , and the green line is the 
corresponding multiscale solution given by formula ( 13.28b . The plots are 
given for different times near the point of gradient catastrophe (x c , t c ) 
of the Hopf solution. Here x c ~ -1.524, t c ~ 0.216. 



which is subject to the asymptotic conditions 1221 

(1.2) A ~ v/-z/2forz -> -oo, 

and 



(1.3) 



A ~ Ai(z) for z — > oo, 



where Ai(z) is the Airy function. Numerically we will consider equation ( II. U on a finite 
interval [zj, z r ] (typically [—10, 10]). The asymptotic solution near ±oo, which will be 
discussed in more detail below, is truncated in a way that the truncation error at zi, z r 
is below 10~ 10 . At these points we impose the values following from the asymptotic 
solutions as boundary conditions, namely 



(1.4) 



A{zi) = y/-Zi/2 
1 

A(z r ) = 



1 



8^2 



{-zi) 



-5/2 



73 



128V2 



-11/2 



2y / 7TZr' /4 



exp 



To solve equation ( II. \\ for z G [z/ , z r ] we use spectral methods since they allow for 
an efficient numerical approximation of high accuracy. We map the interval [zi , z r ] with a 
linear transformation z — > x to the interval 7 = [—1,1] and expand A there in Chebychev 
polynomials. 
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Let us briefly summarize the Chebychev approach, for details see e.g. J3] [161 [T7) ■ The 
Chebyshev polynomials T n (x) are defined on the interval / by the relation 

T n (cos(t)) = cos(ni) , where x = cos(t) , t G [0, tt] . 

A function / on / is approximated via Chebychev polynomials, / w S^Lo a ™^"( x ) 
where the spectral coefficients o„ are obtained by the conditions f(x{) = Xm=o a nT n (xi), 
I = 0, . . . , AT. This approach is called a collocation method. If the collocation points are 
chosen to be xi = cos(irl/N), the spectral coefficients follow from / via a Discrete Cosine 
Transform (DCT) for which fast algorithms exist. We use here a DCT within Matlab. A 
recursive relation for the derivative of Chebychev polynomials implies that the action of 
the differential operator d x on f(x) leads to an action of a matrix D on the vector of 
the spectral coefficients a n . Thus we express A(x) in terms of Chebychev polynomials, 
A(x) — J2n=o A n T n (x) (we typically work with N = 128), and the coefficients of d x A 
in terms of Chebychev polynomials are determined then via DA. 

To solve equation ( II . lb on the interval [zi, z r ], we use an iterative approach, 

(1.5) A 1l+hzz = zA n + 2Al, neN. 

We start with A\{z) = (l + z 2 ) 1//4 /(l + cxp(z))/v / 2- In each step of the iteration we solve 
equation ( 11.51 ) for A n+ i with the boundary conditions ( 11.41 ). The boundary conditions are 
imposed with a r-method: the last two rows of the matrix D 2 for the second derivative are 
replaced with the boundary conditions at a; = ±1. Since T„(±l) = (±1)", the resulting 
matrix L which will be inverted in each step of the iteration, has only 1 and —1 in the last 
two rows and is thus better conditioned than the matrix D 2 . It turns out that the iteration is 
unstable if no relaxation is used. We thus define A n+ i = [iL~ l (zA n + 2Af l ) + (1 — (i)A n 
with fi = 0.009. With this choice of the parameters, the iteration converges. It is stopped 
when the difference between A n+ i and A n is of the order of machine precision (Matlab 
works internally with a precision of the order of 10 -16 ; due to rounding errors machine 
precision is typically limited to the order of 10" 14 ). The solution is shown in Fig.fTTI 
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Figure 1 1 . Hastings-McLeod solution of the Painleve II equation. 

To test the accuracy of the solution we plot in Fig. [12] the quantity PuA as computed 
with spectral methods on the collocation points. It can be seen that the error is biggest 
on the boundary which is even more obvious from Fig. [13] The found solution is also 
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compared to a numerical solution with a standard ode solver as bvp4c in Matlab. The 
solutions agree within the limits of numerical precision. 
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Figure 12. Accuracy of the solution of the Painleve II equation. 
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FIGURE 1 3 . Accuracy of the solution of the Painleve II equation near 
the boundary of the considered interval. 

For general values of z the solution is obtained as follows: for values of z G [zi, z r ] 
they follow from the spectral data via A(z) = J2n=o A n T n (z). Notice that the accuracy 
of the solution is best on the collocation points, but we can expect it to be of the order of 
at least 10 -6 even at points z in between. For values of z < zi, we use the approximation 
A(z) = y/-z/2 - (-z)" 5 / 2 /8/V2 - z^ 11 / 2 , for values of z > z r , we use the 

approximation A(z) = cxp (— |z 3//2 ) / {2^/ttz 1 ^). This provides a global approximation 
to the solution with an accuracy of the order of 1CP 6 and better, which is sufficient for our 
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purposes. Higher precision can be reached within the used approach without problems: 
one can either increase the values of — z; and z r and use a higher number of polynomials, 
or use higher order terms in the asymptotic solution of A for z — ► ±00. 
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